Resonance fluorescence localization microscopy with sub-wavelength resolution

ABSTRACT

The resonance fluorescence spectrum of a number of two-level atoms is driven by a gradient coherent laser field. In the weak dipole-dipole interaction region (separation less than λ/50), a very strong laser field may be applied such that the Rabi frequency is much larger than the dipole-dipole interaction energy. From the spectrum, the positions of each atom may be determined by just a few measurements. This sub-wavelength microscopy scheme is entirely based on far-field technique and it does not require point-by-point scanning, which makes the method more time-efficient. When two atoms are very close to each other (less than λ/50), the position information for each atom may still be obtained with very high accuracy provided that they are not too close to other atoms. The method may be extended to an arbitrarily large region without requiring more peak laser power and only a few measurements are required.

CROSS-REFERENCE TO RELATED APPLICATIONS

None

STATEMENT REGARDING FEDERALLY SPONSORED RESEARCH OR DEVELOPMENT

Not Applicable

BACKGROUND OF THE INVENTION

1. Field of the Invention

This invention relates to optical microscopy. More particularly, it relates to far-field resonance fluorescence localization microscopy.

2. Description of the Related Art Including Information Disclosed Under 37 CFR 1.97 and 1.98

The resolution limit of traditional far-field optical microscopy is about half the wavelength of the light [see, e.g., E. Abbe, Arch. Mikr. Anat. 9, 413 (1873) and L. Rayleigh, Philos. Mag. 47, 81 (1874)]. To achieve better resolution, one must switch to shorter wavelengths (e.g., electron microscopy) which is usually invasive to the system [see, e.g., A. Diaspro (ed), Nanoscopy and Multidimensional Optical Fluorescence Microscopy (CRC Press, Boca Raton, 2010)]. Near-field scanning microscopy can obtain optical imaging with sub-diffraction resolution [see, e.g., G. Binnig, C. F. Quate, and C. Gerber, Phys. Rev. Lett. 56, 930 (1986) and C. Hettich et al., Science 298, 385 (2002)], but due to the surface bound nature it is limited in application. Two-photon fluorescence microscopy first was developed to achieve a higher resolution than classical one-photon fluorescence microscopy in the far field [see, e.g., W. Denk, J. H. Strickler, and W. W. Webb, Science 248, 73 (1990) and J. H. Strickler and W. W. Webb, Proc. SPIE 1398, 107 (1991)]. Stimulated emission depletion (STED) and the related concept of ground-state depletion microscopy are then developed to overcome the far-field diffraction limit in fluorescence microscopy [see, e.g., S. W. Hell and J. Wichmann, Opt. Lett. 19, 780 (1994) and S. W. Hell and M. Kroug, Appl. Phys. B 60, 495 (1995)]. Space-dependent dark states are also proposed to achieve subwavelength resolution [see, e.g., G. S. Agarwal and K. T. Kapale, J. Phys. B 39, 3437 (2006) and S. Bretschneider, C. Eggeling, and S. W. Hell, Phys. Rev. Lett. 98, 218103 (2007)]. However, realization of these schemes is based on point-by-point scanning and is time consuming. Coherent Rabi oscillations may also be employed to break the diffraction limit [see, e.g., Z. Liao, M. Al-Amri, and M. S. Zubairy, Phys. Rev. Lett. 105, 183601 (2010) and C. Shin et al., J. Lumin. 130, 1635 (2010)], but the effect of dipole-dipole interaction has not been well discussed. Another method based on resonance fluorescence is able to measure the separation of two interacting atoms with subwavelength resolution [see, e.g., J-T Chang, J. Evers, M. O. Scully, and M. S. Zubairy, Phys. Rev. A 73, 031803(R) (2006) and Q. Sun, M. Al-Amri, M. O. Scully, and M. S. Zubairy, Phys. Rev. A 83, 063818 (2011)].

A question remains whether the locations of multiple atoms can be determined with sub-wavelength resolution even when dipole-dipole interaction is involved.

BRIEF SUMMARY OF THE INVENTION

An evaluation of the resonance fluorescence spectrum of a number of two-level atoms driven by a gradient coherent laser field shows that the positions of atoms can be determined from the spectrum even when the atoms locate within the sub-wavelength range and the dipole-dipole interaction is significant. This far-field resonance fluorescence localization microscopy method does not require point-by-point scanning and is therefore more time-efficient. The present invention also includes a method to extract the position information in an extended region without requiring additional peak laser power. One particular embodiment of the invention is a method that permits 2D imaging.

BRIEF DESCRIPTION OF THE SEVERAL VIEWS OF THE DRAWING(S)

FIG. 1 depicts a scheme for resonance fluorescence microscopy wherein x_(i) is the position of the ith atom.

FIG. 2 is a dressed state picture for three interacting atoms. |±±±n

=(|bn₁

±|an₁−1

)

(|bn₂

±|an₂−1

)

(|bn₃

±|an₃−1)) where n₁+n₂+n₃=n.

FIG. 3 depicts resonance fluorescence spectra for: a) x₁=0.3λ, x₂=0.5λ, x₃=0.7λ and Ω(x)=100γx/λ; b) x₁=0.45λ, x₂=0.5λ, x₃=0.56λ and Ω(x)=1500γx/λ; c) x₁=0.485λ, x₂=0.5λ, x₃=0.6λ and Ω(x)=1γx/λ; and, d) same as c but with Ω(x)=200γx/λ.

FIG. 4 illustrates a schematic setup for imaging atoms in an extended region based on resonance fluorescence localization microscopy.

FIG. 5 depicts an example of 2D resonance fluorescence localization microscopy according to the invention.

DETAILED DESCRIPTION OF THE INVENTION

The present invention utilizes collective resonance fluorescence to provide the spatial information of a multi-atom system. With this system, far-field resonance fluorescence localization microscopy (RFLM) can be performed with sub-wavelength resolution.

General Formalism for N-Atom Resonance Fluorescence Localization Microscopy

For simplicity, first consider identical atoms located in a line along the x axis. Such a configuration is illustrated in FIG. 1. Two, strong, linearly-polarized laser fields with wavelength λ shine on these atoms from opposite directions and form a standing wave 10. Assume that the polarization orientation is in ŷ direction and the frequency is resonant with the two-level atoms. Also assume that the atoms do not move and they locate within one wavelength. This assumption is valid for the following types of situations: atoms that are trapped by an optical lattice, quantum dots, NV-centers in diamond and the like. The standing wave may be stretched where the sample is located within the approximately linear region between node and anti-node [see, e.g., Z. Liao, M. Al-Amri, and M. S. Zubairy, Phys. Rev. Lett. 105, 183601 (2010) and A. V. Gorshkov, L. Jiang, M. Greiner, P. Zoller, and M. D. Lukin, Phys. Rev. Lett. 100, 093005 (2008)]. In this region, E(x)=E₀x/λ. The resonance fluorescence photons 12 emitted by the system may be monitored with a detector in the {circumflex over (z)} direction. The resonance fluorescence spectrum encodes the spatial information of the systems from which the position of each atom can be determined.

The Hamiltonian of the system and the field is: H=H _(A) +H _(F) +H _(AF) +H _(dd),  (1) where H_(A)=ℏω₀Σ_(i=1) ^(N) S_(i) ^(z) is the energy of the atoms, with ω₀ being the level separation and S_(i) ^(z) is the z component of the spin operator [see, e.g., J-T Chang, J. Evers, M. O. Scully, and M. S. Zubairy, Phys. Rev. A 73, 031803(R) (2006) and T. G. Rudolph, Z. Ficek, and B. J. Dalton, Phys. Rev. A 52, 636 (1995)]. H_(F)=ℏω₀a^(†)a is the total energy of the photons, where a(a^(†)) is the annihilation (creation) operator of the photon; H_(AF)=(ℏ/2)Σ_(i=1) ^(N) g_(i)(S_(i) ⁺a+S_(i) ⁻a^(†)) is the interaction between the atoms and the field, with S_(i) ⁺(S_(i) ⁻) being the raising (lowering) operator on the ith atom, and coupling constant g_(i)=gx_(i)/λ and g=μ(2ω₀/ℏε₀V)^(1/2) (μ is the transition dipole moment between ground state and excited state); H_(dd)=ℏΣ_(i≠j)Ω_(ij)(S_(i) ⁺S_(j) ⁻+S_(i) ⁻S_(j) ⁺) is the dipole-dipole interaction energy. All transition dipole moments are polarized in the y direction and the dipole-dipole interaction energy Ω_(ij) is given by:

$\begin{matrix} {{\Omega_{ij} = {\frac{3\gamma}{4}\left\lbrack {{- \frac{\cos\left( {kx}_{ij} \right)}{{kx}_{ij}}} + \frac{\sin\left( {kx}_{ij} \right)}{\left( {kx}_{ij} \right)^{2}} + \frac{\cos\left( {kx}_{iij} \right)}{\left( {kx}_{ij} \right)^{3}}} \right\rbrack}},} & (2) \end{matrix}$ with 2γ=4ω₀ ³d₀ ²/(3ℏc³) being the single-atom spontaneous decay rate, k=ω₀/c (c is speed of light) and x_(ij) is the distance between atoms [see, e.g., Z. Ficek and S. Swain, Quantum Interference and Quantum Coherence: Theory and Experiment (Springer, N.Y., 2004)]. The Rabi frequency for the ith atom is given by Ω_(i)=g_(i)√{square root over (n)} (or μE₀x_(i)/ℏλ) where n is the photon number.

For Ω_(i)>>Ω_(ij)

If the dipole-dipole interaction energy is not very strong, a strong laser field may be applied such that Ω_(i)>>Ω_(ij). In this case, the collective resonance fluorescence spectrum of a multi-atom system may be analytically evaluated based on a dressed state picture [see, e.g., C. Cohen-Tannoudji and S. Reynaud, J. Phys. B: Atom. Molec. Phys. 10, 345 (1977) and H. S. Freedhoff, Phys. Rev. A 19, 1132 (1979)]. Let H₀=H_(A)+H_(F)+H_(AF) and treat H_(dd) as a perturbation term. The eigenvalues and eigenfunctions of H₀ are given by:

$\begin{matrix} {\mspace{79mu}{{E_{\alpha,n}^{0} = {{\left( {n - \frac{N}{2}} \right){\hslash\omega}_{0}} + {\frac{\hslash}{2}{\sum\limits_{i = 1}^{N}{\chi_{i}^{\alpha}\Omega}}}}},{\left. {\alpha,n} \right\rangle = {\frac{1}{\sqrt{2^{N}}}\left( {\left. {b^{N},n} \right) + {\sum\limits_{i = 1}^{N}{\chi_{i}^{\alpha}\left. {{a_{i}b^{N - 1}},{n - 1}} \right\rangle}} + {\sum\limits_{i \neq j}{\chi_{i}^{\alpha}\chi_{j}^{\alpha}\left. {{a_{i}a_{j}b^{N - 2}},{n - 2}} \right\rangle}} + \ldots + {\prod\limits_{i = 1}^{N}\;{\chi_{i}^{\alpha}\left. {a^{N},{n - N}} \right\rangle}}} \right)}},}} & (3) \end{matrix}$ where N is the number of atoms, α=1, 2, . . . , 2^(N), |a_(i)b^(N-1)

means that ith atom is in the excited state |a

while other N−1 atoms are in the ground state |b

, and χ_(i) ^(α) is a constant which can be either +1 or −1. Counting the dipole-dipole interaction term as a perturbation, the eigenenergy is shifted by

$\begin{matrix} {\Delta_{\alpha,n} = {\frac{\hslash}{2}{\sum\limits_{i \neq j}{\chi_{i}^{\alpha}\chi_{j}^{\alpha}\Omega_{ij}}}}} & (4) \end{matrix}$ and the correction to zeroth order eigenfunction is of the order of Ω_(ij)/Ω_(i) which can be neglected. The sublevel energy E_(α,n)=E_(α,n) ⁰+Δ_(α,n). The pictorial energy level for the dressed state picture is shown in FIG. 2. Coupling of the dressed states to the vacuum results in the system's cascade down the ladder from the α state of one multiplet to the β state of the adjacent multiplet [see, e.g., C. Cohen-Tannoudji and S. Reynaud, J. Phys. B: Atom. Molec. Phys. 10, 345 (1977)]. The corresponding transition frequency is ω_(αβ)=(E_(α,n)−E_(β,n-1))/ℏ, where α, β=1, 2, . . . , 2^(N).

The spectrum of resonance fluorescence can be evaluated by:

$\begin{matrix} {{{S(\omega)} \propto {{Re}\left\lbrack {\int_{0}^{\infty}{{\mathbb{d}{\tau\mathbb{e}}^{\mathbb{i}\omega\tau}}{\lim\limits_{t\rightarrow\infty}\left\langle {{D^{+}(t)}{D^{-}\left( {t + \tau} \right)}} \right\rangle}}} \right\rbrack}},} & (5) \end{matrix}$ where D⁺ and D⁻ are the raising and lowering parts of the total atomic dipole operator [see, e.g., C. Cohen-Tannoudji and S. Reynaud, J. Phys. B: Atom. Molec. Phys. 10, 345 (1977) and H. S. Freedhoff, Phys. Rev. A 19, 1132 (1979)]. The lowering part may be written as D⁻=Σ_(αβn) d_(αβ) ⁻|β, n−1

α, n|=Σ_(αβ)D_(αβ) ⁻, where d_(αβ) ⁻ is the dipole matrix element of the transition from |α, n

to |β, n−1

and it is defined by d_(αβ) ⁻=

β, n−1|Σ_(i=1) ^(N)S_(i) ⁻|α, n

and D_(αβ) ⁻=Σ_(n)d_(αβ) ⁻|β, n−1

α, n|.

The two-time correlation function in Eq. (6) may be written as:

D⁺(t)D ⁻(t+τ)

=Σ_(α≠β)

D⁺(t)D _(αβ) ⁻(t+τ)+Σ_(α)

D⁺(t)D _(αα) ⁻(t+τ)  (6) where the first term corresponds to the sideband spectrum while the second term corresponds to the central peak. According to the quantum regression theorem [see, e.g., M. O. Scully and M. S. Zubairy, Quantum Optics (Cambridge University Press, Cambrige, 1997)], the two-time correlation function

D⁺(t)D_(αβ) ⁻(t+τ)

satisfies the same equation of motion as the single time average

D_(αβ) ⁻(t)

. The dynamics of

D_(αβ) ⁻(t)

can be calculated from the master equation:

$\begin{matrix} {\frac{\mathbb{d}\left\langle {D_{\alpha\beta}^{-}(t)} \right\rangle}{\mathbb{d}t} = {{d_{\alpha\beta}^{-}\frac{\mathbb{d}\rho_{\alpha\beta}^{-}}{\mathbb{d}t}} = {d_{\alpha\beta}^{-}\left\{ {{\frac{- i}{\hslash}\left\lbrack {H,\rho} \right\rbrack}_{\alpha\beta} - \left( {L\;\rho} \right)_{\alpha\beta}} \right\}}}} & (7) \end{matrix}$ where ρ_(αβ) ⁻=

α, n|ρ|β, n−1

, L=Σ_(i,j=1) ^(N)γ_(ij)(S_(i) ⁺S_(j) ⁻ρ+ρS_(i) ⁺S_(j) ⁻−2S_(j) ⁻ρS_(i) ⁺) is the relaxation operator with γ_(ii) being the decay rate of atom i and γ_(ij) being the cross damping rate. For the sidebands, one may expand (Lρ)_(αβ)=Γ_(αβ)ρ_(αβ)+ . . . and, from Eq. (8), one obtains:

$\begin{matrix} {{{\frac{\mathbb{d}}{\mathbb{d}t}\left\langle {D_{\alpha\beta}^{-}(t)} \right\rangle} \simeq {\left( {{\mathbb{i}\omega}_{\alpha\beta} - \Gamma_{\alpha\beta}} \right)\left\langle {D_{\alpha\beta}^{-}(t)} \right\rangle}},} & (8) \end{matrix}$ where the non-resonance terms on the right hand side in the secular approximation have been neglected. For the central peak, as all |α, n

→|α, n−1

, α=1, . . . , 2^(N) have the same transition frequency, they couple to each other and one may expand (Lρ)_(αα)=Σ_(β)Γ_(αβ)′ρ_(ββ)+ . . . . From Eq. (8) one obtains:

$\begin{matrix} {{\frac{\mathbb{d}}{\mathbb{d}t}\left\langle {D_{\alpha\alpha}^{-}(t)} \right\rangle} = {{{\mathbb{i}\omega}_{0}\left\langle {D_{\alpha\alpha}^{-}(t)} \right\rangle} - {d_{\alpha\alpha}^{-}{\sum\limits_{\beta}\;{\Gamma_{\alpha\beta}^{\prime}{\frac{\left\langle {D_{\beta\beta}^{-}(t)} \right\rangle}{d_{\beta\beta}^{-}}.}}}}}} & (9) \end{matrix}$

According to the quantum regression theorem and Eq. (6), the spectrum is given by:

$\begin{matrix} {{{S\left( {\overset{\rightarrow}{R},\omega} \right)} = {{{S^{0}(\omega)} + {S^{\pm}(\omega)}} \propto {{{Re}\left\lbrack {\int_{0}^{\infty}{{\mathbb{d}{\tau\mathbb{e}}^{{{\mathbb{i}}{({\omega - \omega_{0}})}}\tau}}{\sum\limits_{\alpha,\beta}{{d_{\alpha\alpha}^{-}\left( {\mathbb{e}}^{{- \Gamma^{\prime}}\tau} \right)}_{\alpha\beta}d_{\beta\beta}^{+}}}}} \right\rbrack} + {\sum\limits_{\alpha \neq \beta}\frac{{d_{\alpha\beta}^{-}}^{2}\Gamma_{\alpha\beta}}{\left( {\omega - \omega_{\alpha\beta}} \right)^{2} + \Gamma_{\alpha\beta}^{2}}}}}},} & (10) \end{matrix}$ where the first term yields the central peak spectrum and the second term gives the sideband spectrum.

For the zeroth order wavefunction (Eq. (4) or equivalently Eq. (21)), the transition dipole is given by:

$\begin{matrix} {d_{\alpha\;\beta}^{-} = {\frac{1}{2^{N}}{\sum\limits_{i = 1}^{N}{\left\{ {\chi_{i}^{\beta}{\prod\limits_{k \neq i}\;\left\lbrack {1 + {\chi_{k}^{\alpha} \cdot \chi_{k}^{\beta}}} \right\rbrack}} \right\}.}}}} & (11) \end{matrix}$

There are three possible cases:

Case 1. β=α, d_(αβ) ⁻=Σ_(i=1) ^(N)χ_(i) ^(α)/2 which contributes to the central peak ω=ω₀;

Case 2. β=α^(p) (α^(p) is a state such that E_(α) ⁰ and E_(α) _(p) ⁰ have different sign only in pth term), d_(αβ) ⁻=χ_(p) ^(α)/2 which contributes to the sidebands ω_(αα) _(p) =ω₀+χ_(p) ^(α)Ω_(p)+Σ_(k≠p)χ_(p) ^(α)χ_(k) ^(α)Ω_(pk);  (12)

From this equation, it is apparent that the positive sideband peaks can be divided into N groups: Ω_(p)+Σ_(k≠p)±Ω_(pk), p=1, . . . , N. Averaging over the frequencies of each group provides the Rabi frequencies Ω_(p) from which the positions of the atoms may be determined. The error is on the order of Ω_(ij) ²/Ω_(i) ²<<1. This is the method of the present invention for optical microscopy. In an experiment, it may not be known which peak belongs to which group. However, if the gradient of laser field is changed [by an amount], the relative Rabi frequencies for different atoms change which causes the separations between different groups of spectrum shift. Because the dipole-dipole interactions do not change, the splitting between peaks belong to the same group will not change. From this phenomena, peaks belonging to different groups can be identified.

Case 3. α≠β and more than one term of E_(α) ⁰ and E_(β) ⁰ have different signs, d_(αβ) ⁻=0 which corresponds to the forbidden transition.

The method described above is valid under the conditions: Ω_(i)>>Ω_(ij) and |Ω_(i)−Ω_(j)|−2Ω_(ij)>>γ. Assuming that γ˜10⁸ Hz and the maximum Rabi frequency is 10¹³ Hz, then the smallest distance that can be resolved using this method is about λ/50.

For Ω_(ij)>>Ω_(i)

When there are two atoms in the sample whose distance is very close (e.g. r_(ij)<λ/50), the condition Ω_(i)>>Ω_(ij) cannot be satisfied. The positions of these two atoms cannot be localized based on the method described above. However, there are ways to extract the position information of the two atoms if they are far away from other atoms (e.g., greater than λ/10). In this case a weak gradient field may be applied such that Ω_(i), Ω_(j)<<Ω_(ij). If Rabi frequency Ω_(i)˜γ, there are only two sideband peaks located at ω₀±Ω_(ij) [see, e.g., J-T Chang, J. Evers, M. O. Scully, and M. S. Zubairy, Phys. Rev. A 73, 031803(R) (2006) and T. G. Rudolph, Z. Ficek, and B. J. Dalton, Phys. Rev. A 52, 636 (1995)]. Therefore, the dipole-dipole interaction energy Ω_(ij) can be determined from the resonance fluorescence spectrum. According to Eq. (2), the distance r_(ij) between these two, close atoms can then be determined.

In one experiment, the gradient field strength was increased to a medium value (for example ˜100γ) which is still much less than dipole-dipole interaction energy. At this point, each sideband peak split into two peaks [see, id.]. For a positive sideband, it splits into two peaks:

$\begin{matrix} {{\omega_{+}^{1} = {\omega_{0} + \Omega_{ij} + \frac{\left. {\Omega_{i}\Omega_{j}} \right)^{2}}{2\Omega_{ij}}}},} & (13) \\ {\omega_{+}^{2} = {\omega_{0} + \Omega_{ij} + {\frac{\Omega_{i}\Omega_{j}}{\Omega_{ij}}.}}} & (14) \end{matrix}$

From each equation and the relationship between Ω_(i) and Ω_(j), a value for the positions r_(i) and r_(j) can be calculated. Because there are two equations, there may be two results. One of the results may be used as the position[s] of the two atoms or the two results may be averaged to obtained the positions of the two atoms.

Linewidth

In general, the linewidth of the emitted radiation is difficult to calculate exactly. However, the linewidth may be evaluated approximately in some cases (discussed infra). From Eq. (2), when the distance between two atoms is about λ/10, the dipole-dipole interaction energy is comparable to the linewidth of the sideband spectrum of the independent atoms. Therefore, λ/10 may be set as a threshold and the linewidth evaluated. First, when all the atoms have distances much larger than λ/10, the dipole-dipole interaction energies are much smaller than the sideband spectrum linewidth of independent atoms and the dipole-dipole interaction energies may be neglected. For independent atoms, the linewidth of the sideband spectrum is 3γ/2 [see, e.g., H. S. Freedhoff, Phys. Rev. A 19, 1132 (1979) and M. O. Scully and M. S. Zubairy, Quantum Optics (Cambridge University Press, Cambrige, 1997)], i.e., Γ_(αβ)=3γ/2 in Eq. (11). Second, when all the atoms have distances much smaller than λ/10, all dipole-dipole interaction energies Ω_(ij) are larger than 3γ/2 and the overlapped sideband spectrum splits. One can calculate (Lρ)_(αβ)≈(N/2+1)γρ_(αβ)+ . . . , from which it may be seen that Γ_(αβ)≈(N/2+1)γ. The spectrum width is about (N/2+1)γ which is similar to superradiance [see, e.g., R. H. Dicke, Phys. Rev. 93, 99 (1954)].

For the general case wherein some atoms have distances larger than λ/10 and some atoms have distances smaller than λ/10, the resonance fluorescence has the same transition frequencies but a different linewidth which is about (N_(eff)+2)γ/2 where N_(eff) is the average number of atoms which couple to each other and its value is between 1 and N.

Three-Atom Example

In the following, the resonance fluorescence spectrum of a three-atom system is solved numerically to demonstrate how this localization microscopy works. The dressed state picture is shown in FIG. 2. From Eq. (3), the eigenvalues of the system are:

$E_{1\; n} = {{\left( {n - \frac{3}{2}} \right)\hslash\;\omega_{0}} + {\frac{\hslash}{2}\left( {\Omega_{1} + \Omega_{2} + \Omega_{3} + \Omega_{12} + \Omega_{13} + \Omega_{23}} \right)}}$ $E_{2\; n} = {{\left( {n - \frac{3}{2}} \right)\hslash\;\omega_{0}} + {\frac{\hslash}{2}\left( {{- \Omega_{1}} + \Omega_{2} + \Omega_{3} - \Omega_{12} - \Omega_{13} + \Omega_{23}} \right)}}$ $E_{3\; n} = {{\left( {n - \frac{3}{2}} \right)\hslash\;\omega_{0}} + {\frac{\hslash}{2}\left( {\Omega_{1} - \Omega_{2} + \Omega_{3} - \Omega_{12} + \Omega_{13} - \Omega_{23}} \right)}}$ $E_{4\; n} = {{\left( {n - \frac{3}{2}} \right)\hslash\;\omega_{0}} + {\frac{\hslash}{2}\left( {\Omega_{1} + \Omega_{2} - \Omega_{3} + \Omega_{12} - \Omega_{13} - \Omega_{23}} \right)}}$ $E_{5\; n} = {{\left( {n - \frac{3}{2}} \right)\hslash\;\omega_{0}} + {\frac{\hslash}{2}\left( {{- \Omega_{1}} - \Omega_{2} + \Omega_{3} + \Omega_{12} - \Omega_{13} - \Omega_{23}} \right)}}$ $E_{6\; n} = {{\left( {n - \frac{3}{2}} \right)\hslash\;\omega_{0}} + {\frac{\hslash}{2}\left( {{- \Omega_{1}} + \Omega_{2} - \Omega_{3} - \Omega_{12} + \Omega_{13} - \Omega_{23}} \right)}}$ $E_{7\; n} = {{\left( {n - \frac{3}{2}} \right)\hslash\;\omega_{0}} + {\frac{\hslash}{2}\left( {\Omega_{1} - \Omega_{2} - \Omega_{3} - \Omega_{12} - \Omega_{13} + \Omega_{23}} \right)}}$ $E_{8\; n} = {{\left( {n - \frac{3}{2}} \right)\hslash\;\omega_{0}} + {\frac{\hslash}{2}\left( {{- \Omega_{1}} - \Omega_{2} - \Omega_{3} + \Omega_{12} + \Omega_{13} + \Omega_{23}} \right)}}$

For weak dipole-dipole interactions, according to Eq. (12) the nonzero transition dipoles are: d₁₁ ^(±)=3/2, d₂₂ ^(±)=d₃₃ ^(±)=d₄₄ ^(±)=½, d₅₅ ^(±)=d₆₆ ^(±)=d₇₇ ^(±)=−½, d₈₈ ^(±)=−3/2; d₁₂ ⁻=d₁₃ ⁻=d₁₄ ⁻=d₂₅ ⁻=d₂₆ ⁻=d₃₅ ⁻=d₃₇ ⁻=d₄₇ ⁻=d₅₈ ⁻=d₆₈ ⁻=d₇₈ ⁻=½; d₂₁ ⁻=d₃₁ ⁻=d₄₁ ⁻=d₅₂ ⁻=d₆₂ ⁻=d₅₃ ⁻=d₇₃ ⁻=d₇₄ ⁻=d₈₅ ⁻=d₈₆ ⁻=d₈₇ ⁻=−½. One may also calculate (Lρ)=_(αα)=(3γ/2)ρ_(αα)−(γ/2)Σ_(α′)ρ_(α′α′) where α→α′ is allowed sideband transition. For example, if α=1, then α′=2,3,4. Thus,

$\Gamma^{\prime} = {\frac{\gamma}{2}\begin{pmatrix} 3 & {- 1} & {- 1} & {- 1} & 0 & 0 & 0 & 0 \\ {- 1} & 3 & 0 & 0 & {- 1} & {- 1} & 0 & 0 \\ {- 1} & 0 & 3 & 0 & {- 1} & 0 & {- 1} & 0 \\ {- 1} & 0 & 0 & 3 & 0 & {- 1} & {- 1} & 0 \\ 0 & {- 1} & {- 1} & 0 & 3 & 0 & 0 & {- 1} \\ 0 & {- 1} & 0 & {- 1} & 0 & 3 & 0 & {- 1} \\ 0 & 0 & {- 1} & {- 1} & 0 & 0 & 3 & {- 1} \\ 0 & 0 & 0 & 0 & {- 1} & {- 1} & {- 1} & 3 \end{pmatrix}}$

From Eq. (9), one may obtain the central peak spectrum:

$\begin{matrix} {{S^{0}\left( {\overset{\rightarrow}{R},\omega} \right)} \propto {{\frac{3\gamma}{2}\left\lbrack {\frac{1}{\left( {\omega - \omega_{0}} \right)^{2} + \gamma^{2}} + \frac{18}{\left( {\omega - \omega_{0}} \right)^{2} + {4\gamma^{2}}} + \frac{9}{\left( {\omega - \omega_{0}} \right)^{2} + {9\gamma^{2}}}} \right\rbrack}.}} & (15) \end{matrix}$

For the sidebands, the eigenenergy and eigenvectors of the dressed system may be numerically solved and the allowed transition frequencies calculated. In an actual experiment, the sideband peaks of the spectrum were measured. According to Eq. (13), there are three groups of allowed sideband spectrum on the positive side: ω₀+Ω₁±Ω₁₂±Ω₁₃, ω₀+Ω₂±Ω₁₂±Ω₂₃ and ω₀+Ω₃±Ω₁₃±Ω₂₃. The Rabi frequencies for each atom may be obtained by averaging over each group of the spectrum. For example, if there are three atoms and their positions are x₁=0.3λ, x₂=0.5λ and x₃=0.7λ, the separation is λ/5. A gradient electric field is applied such that Ω(x)=100γx/λ. The resonance fluorescence spectrum is shown in FIG. 3 a. From the spectrum, the sideband frequencies are Ω₁=(30.00±1.82)γ, Ω₂=(50.00±1.82)γ, Ω₃=(70.00±1.72)γ. The positions of the atoms may be determined as: x₁=(0.300±0.002)λ, x₂=(0.500±0.002)λ, and x₃=(0.700±0.002)λ which match the previously set parameters very well.

As a second example, consider that the three atoms are located at positions x₁=0.45λ, x₂=0.5λ and x₃=0.56λ. The shortest distance is λ/20. A strong gradient electric field is applied such that Ω(x)=1500γx/λ. The resonance fluorescence spectrum is shown in FIG. 3 b. From the spectrum, the sideband frequencies may be determined as shown in the following table (γ):

648.4 ± 2.7 652.4 ± 2.7 694.3 ± 2.8 698.3 ± 2.8 715.5 ± 2.6 741.3 ± 2.5 761.5 ± 2.7 787.2 ± 2.5 825.7 ± 2.7 829.7 ± 2.7 851.4 ± 2.8 855.3 ± 2.8

Then, the gradient laser field was slightly increased such that Ω(x)=1700γx/λ and the spectrum peaks are shown in the following table (γ):

693.6 ± 2.7 697.4 ± 2.7 739.5 ± 2.8 743.3 ± 2.8 765.4 ± 2.6 791.2 ± 2.5 811.4 ± 2.7 837.2 ± 2.5 881.7 ± 2.8 885.5 ± 2.8 907.4 ± 2.8 911.3 ± 2.8

Comparing these two tables, it is apparent that the separations between peaks in each column do not change significantly. However, the separations between peaks in different columns change significantly. From this result it can be determined that the spectrum from each column belongs to the same group. Averaging over each column of the first table one may obtain Ω₁=(673.35±2.75)γ, Ω₂=(751.38±2.58)γ, Ω₃=(840.53±2.75)γ. The positions of the atoms may then be determined: x₁=(0.449±0.002)λ, x₂=(0.501±0.002)λ, and x₃=(0.560±0.002)λ which also match the actual positions of the atoms quite well.

If there are two atoms whose distance is smaller than this limit, the method described above may be used. For example, there are three atoms and their positions are x₁=0.485λ, x₂=0.5λ and x₃=0.6λ. The distance between the first atom and the second atom is 0.015λ which is less than λ/50. From Eq. (2), one may calculate that the dipole-dipole interaction energy is 891.92γ which is very large. In this situation, a weak gradient laser field may be applied first such that the corresponding Rabi frequency is Ω(x)=1γx/λ. The resonance fluorescence spectrum may be solved numerically and the result is shown in FIG. 3 c. It can be seen that only two sideband peaks appear and their positions are ω=ω_(o)±891.90γ which match very well with the calculation value. From Eq. (2), one may deduce that the distance between these two, close atoms is Δx₁₂=0.015λ. Then, the gradient laser field may be increased to a medium value, for example Ω(x)=200γx/λ. The numerical result of the corresponding resonance fluorescence spectrum is shown in FIG. 3 d. From the spectrum it can be seen that each sideband peak split into two peaks. For example, the positive sideband peak splits into ω₀+902.72γ and ω₀+913.36γ. According to Eq. (14) and Eq. (15),

$\begin{matrix} {{\Omega_{12} + \frac{\left( {\Omega_{1} + \Omega_{2}} \right)^{2}}{2\Omega_{12}}} = {913.36\gamma}} & (16) \\ {{\Omega_{12} + \frac{\Omega_{1}\Omega_{2}}{\Omega_{12}}} = {902.72\gamma}} & (17) \end{matrix}$ where Ω₁₂=1794.84γ. Assuming that Ω₁=200γx₁/λ and Ω₂=200γ(X₁+Δx₁₂)/λ, one may get x₁=0.4816λ from Eq. (18) and x₁′=0.4837λ from Eq. (19). Averaging these two results produces x ₁=0.483λ which is very close to the value 0.485λ previously set, and the error is about 0.4%. Additionally, there are two sideband peaks observed near the central peaks which read (±120±2.6)γ. They are the resonance fluorescence from the third atom and its position may be determined to be (0.600±0.013)λ which also matches the correct [preselected] value. Extension to Larger Area and Higher Dimensions

Above, it was disclosed how to resolve atoms located within one wavelength. For a region larger than one wavelength, one simple method is to stretch the standing wave with larger periods to cover the entire region. This method is easy to perform, but one disadvantage is that the field intensity increases as the working region increases. If the region is too large, the field will be incredible large. Therefore, for a working region beyond several wavelengths extension, a new way may be needed. One way to extend this limitation is via the divide-and-conquer method. The scheme is shown in FIG. 4. The sample is first exposed to a standing wave denoted by curve 42. The regions designated 44 on the object plane locate in an approximately linear field region while regions 46 do not. The resonance fluorescences are collected by a lens. The fluorescence emitted by regions 44 are focused to detector pixels 48 on the imaging plane, while the fluorescence emitted by regions 46 are focused to detector pixels 49. At this step, only the spectrum of the fluorescence collected by detector pixels 48 are analyzed and the positions of the atoms in regions 44 on the object plane may be determined by the method disclosed above. Then, the standing wave is shifted by a phase λ/2. At this time, regions 46 locate in an approximately linear field while regions 44 do not. Applying a similar process, the positions of the atoms in regions 46 may be determined. Because the image of a point in the object plane is not a point but a small disk which is usually described by the point spread function of the lens, there may be a gap 47 between neighboring detector pixels (48, 49) to ensure that the fluorescence from regions 44 do not shine on detector pixels 49 and vice versa. If there is an optical detecting array for each working region on the order of several wavelengths, this method is possible.

The method of the present invention may also be applied to a 2D image. Such a scheme is shown in FIG. 5. Three steps are needed to obtain the two dimensional spatial information. In the first two steps, a gradient laser field is shined along the x(y) direction, from the resonance fluorescence spectrum a discrete set of x(y) position information of the atoms may be obtained. After that, the positions of the atoms still cannot be determined because all combinations of x values and y values are possible. One should therefore shine a third gradient field from a direction which avoids any two pairs of (x, y). From the third resonance fluorescence spectrum the positions of the atoms may be determined.

The resonance fluorescence spectrum of a number of two-level atoms is driven by a gradient coherent laser field. In the weak dipole-dipole interaction region (separation less than λ/50), a very strong laser field may be applied such that the Rabi frequency is much larger than the dipole-dipole interaction energy. From the spectrum the positions of each atom may be determined by just a few measurements. This subwavelength microscopy scheme is entirely based on far-field technique and it does not require point-by-point scanning, which makes the method more time-efficient. When two atoms are very close to each other (less than λ/50), the position information for each atom may still be obtained with very high accuracy provided that they are not too close to other atoms. The method may be extended to an arbitrary large region without requiring more peak laser power and only a few measurements are required.

In the method of the present invention, uncertainty may be due to the linewidth of the spectrum and the calibration of the light intensity. Lower density of atoms within λ/10 provides narrower linewidth and therefore smaller uncertainty. Good calibration of the light intensity is also required in an experiment in order to extract more precise position information. There is a limitation on the number of atoms within one wavelength which is about 50 in 1D and 2500 in 2D space. Another limitation of the method is that it is still not clear how to extract the spatial information of atoms when more than two atoms are very close to each other (less than λ/50).

Eigenvalues and Eigenvectors

The eigenvalue and eigenvector for H₀ ^(i) are:

$\begin{matrix} {{{E_{i}^{+} = {{\left( {n_{i} - \frac{1}{2}} \right)\hslash\;\omega_{0}} + \frac{{\hslash\Omega}_{i}}{2}}},{\left. \Phi_{i}^{+} \right\rangle = \frac{\left. {g_{i},n_{i}} \right\rangle + \left. {e_{i},{n_{i} - 1}} \right\rangle}{\sqrt{2}}}}{{E_{i}^{-} = {{\left( {n_{i} - \frac{1}{2}} \right)\hslash\;\omega_{0}} - \frac{{\hslash\Omega}_{i}}{2}}},{\left. \Phi_{i}^{-} \right\rangle = \frac{\left. {g_{i},n_{i}} \right\rangle - \left. {e_{i},{n_{i} - 1}} \right\rangle}{\sqrt{2}}}}} & (18) \end{matrix}$ where Ω_(i)=g√{square root over (n_(i))} and n_(i) is the mean number of photons that interact with the ith atom. The eigenvalues of H₀ are just the summation of the eigenvalues of each atom:

$\begin{matrix} {E_{\alpha}^{0} = {{\left( {n - \frac{N}{2}} \right){\hslash\omega}_{0}} + {\frac{\hslash}{2}{\sum\limits_{i}{\chi_{i}^{\alpha}\Omega_{i}}}}}} & (19) \end{matrix}$ where α=1, 2, . . . , 2^(N), n=Σ_(i) n_(i) and χ_(i) ^(α)=±1. The corresponding eigenvectors are:

$\begin{matrix} \left. \left. {\left. \alpha \right\rangle = {{{\frac{1}{\sqrt{2^{N}}}\overset{N}{\underset{i = 1}{\otimes}}\left( {\left. {g_{i},n_{i}} \right\rangle + \chi_{i}^{\alpha}} \right.}e_{i}n_{i}} - 1}} \right\rangle \right) & (20) \end{matrix}$ which is equivalent to Eq. (2).

The perturbation energy due to the dipole-dipole interaction is:

$\begin{matrix} \begin{matrix} {\Delta_{\alpha} = \left\langle {\alpha{H_{dd}}\alpha} \right\rangle} \\ {= {\hslash{\sum\limits_{i \neq j}{\Omega_{ij}\left\langle {\alpha{{{S_{i}^{+}S_{j}^{-}} + {S_{j}^{+}S_{i}^{-}}}}\alpha} \right\rangle}}}} \\ {= {\frac{\hslash}{4}{\sum\limits_{i \neq j}{\Omega_{ij}\left( \left\langle {g_{i}{{+ \left\langle {e_{i}\left. \chi_{i}^{\alpha} \right)\left( \left\langle {g_{j}{{{+ e_{j}}\left. \chi_{j}^{\alpha} \right)}}} \right. \right.} \right.}}} \right. \right.}}}} \\ {\left( {{S_{i}^{+}S_{j}^{-}} + {S_{j}^{+}S_{i}^{-}}} \right)\left( {\left. g_{i} \right\rangle + {\chi_{i}^{\alpha}\left. e_{i} \right\rangle}} \right)\left( {\left. g_{i} \right\rangle + {\chi_{j}^{\alpha}\left. e_{j} \right\rangle}} \right)} \\ {= {\frac{\hslash}{4}{\sum\limits_{i \neq j}{\Omega_{ij}\left( \left\langle {g_{i}{{+ \left\langle {e_{i}\left. \chi_{i}^{\alpha} \right)\left( \left\langle {g_{j}{{+ \left\langle {e_{j}\left. \chi_{j}^{\alpha} \right)} \right.}}} \right. \right.} \right.}}} \right. \right.}}}} \\ {\left( {{\chi_{j}^{\alpha}\left. e_{i} \right\rangle\left. g_{i} \right\rangle} + {\chi_{i}^{\alpha}\left. g_{i} \right\rangle\left. e_{j} \right\rangle}} \right)} \end{matrix} & (21) \end{matrix}$ wherein the photon part has been ignored because the dipole-dipole Hamiltonian only depends on atomic operators.

Linewidth

Assume that |α

→|β

is an allowed transition, i.e., α=β or they differ from each other by only one term. The relaxation term for the |α

→|β

transition is given by: (Lρ)_(αβ)=Σ_(i,j=1) ^(N)γ_(ij)(

α|S_(i) ⁺ S _(j) ⁻Σ|β

+

α|ΣS_(i) ⁺ S _(j) ⁻|β

−2

α|S_(j) ⁻ ρS _(i) ⁺|β

)

For i=j,

$\begin{matrix} \begin{matrix} {{S_{i}^{+}S_{i}^{-}\left. \alpha \right\rangle} = {{\frac{1}{\sqrt{2^{N}}}\underset{k \neq i}{\otimes}\left( {\left. g_{k} \right\rangle + {\chi_{k}^{\alpha}\left. e_{k} \right\rangle}} \right)}\left( {\chi_{i}^{\alpha}\left. e_{i} \right\rangle} \right)}} \\ {= {\frac{1}{2}\left( {\left. \alpha \right\rangle - \left. \alpha^{i} \right\rangle} \right)}} \end{matrix} & (23) \\ \begin{matrix} {{S_{i}^{+}S_{i}^{-}\left. \beta \right\rangle} = {{\frac{1}{\sqrt{2^{N}}}\underset{k \neq i}{\otimes}\left( {\left. g_{k} \right\rangle + {\chi_{k}^{\beta}\left. e_{k} \right\rangle}} \right)}\left( {\chi_{i}^{\beta}\left. e_{i} \right\rangle} \right)}} \\ {= {\frac{1}{2}\left( {\left. \beta \right\rangle - \left. \beta^{i} \right\rangle} \right)}} \end{matrix} & (24) \\ \begin{matrix} \left. {{S_{i}^{+}\left. \alpha \right\rangle} = {{\frac{1}{\sqrt{2^{N}}}\underset{k \neq i}{\otimes}\left( {\left. g_{k} \right\rangle + {\chi_{k}^{\alpha}\left. e_{k} \right\rangle}} \right)}\left. e_{i} \right\rangle}} \right) \\ {= {\chi_{i}^{\alpha}}^{\frac{{\alpha\rangle} - {\alpha^{i}\rangle}}{2}}} \end{matrix} & (25) \\ \begin{matrix} \left. {{S_{i}^{+}\left. \beta \right\rangle} = {{\frac{1}{\sqrt{2^{N}}}\underset{k \neq i}{\otimes}\left( {\left. g_{k} \right\rangle + {\chi_{k}^{\beta}\left. e_{k} \right\rangle}} \right)}\left. e_{i} \right\rangle}} \right) \\ {= {\chi_{i}^{\beta}}^{\frac{{\beta\rangle} - {\beta^{i}\rangle}}{2}}} \end{matrix} & (26) \end{matrix}$ where |α^(i)

and |α

have different sign only on the ith term and |β^(i)

and |β

have different sign on the ith term. There are N−1 pairs which satisfy this condition. From (26-29), one may obtain:

$\begin{matrix} {\sum\limits_{i}{\gamma_{ii}\left( \left\langle {\alpha{{{S_{i}^{+}S_{i}^{-}\rho\left. \beta \right\rangle} + \left\langle {{\alpha\left. {{\rho\; S_{i}^{+}S_{i}^{-}\left. \beta \right\rangle} - {2\left\langle {\alpha{{S_{i}^{-}\rho\; S_{i}^{+}}}\beta} \right\rangle}} \right)} = {\frac{\gamma}{2}{\sum\limits_{i}\left\{ {{\left\lbrack {2 - {\chi_{i}^{\alpha}\chi_{i}^{\beta}}} \right\rbrack\rho_{\alpha\beta}} - {\left\lbrack {1 - {\chi_{i}^{\alpha}\chi_{i}^{\beta}}} \right\rbrack\rho_{\alpha^{i}\beta}} - {\left\lbrack {1 - {\chi_{i}^{\alpha}\chi_{i}^{\beta}}} \right\rbrack\rho_{{\alpha\beta}^{i}}} - {\chi_{i}^{\alpha}\chi_{i}^{\beta}\rho_{\alpha^{i}\beta^{i}}}} \right\}}}} \right.}}} \right. \right.}} & (27) \\ {\mspace{79mu}{= \left( \begin{matrix} {{\frac{N_{\gamma}}{2}\rho_{\alpha\beta}} - {\frac{\gamma}{2}{\sum\limits_{i}\rho_{\alpha^{i}\beta^{i}}}}} \\ {{\frac{\left( {N + 2} \right)\gamma}{2}\rho_{\alpha\beta}} - {\frac{\gamma}{2}{\sum\limits_{i}\rho_{\alpha^{i}\beta^{i}}}} - \rho_{\beta\alpha} - {2\rho_{\beta\beta}} - {2\rho_{\alpha\alpha}}} \end{matrix} \right.}} & (28) \end{matrix}$ in which the first equation in Eq. (28) is for α=β, while the second equation is for α≠β but only differs in one term. The last three terms in the second equation may be ignored in the secular approximation because they have a different transition frequency from ω_(αβ). The survival of the second term depends on the coupling of the system. If all dipole-dipole interaction energies are small and can be neglected, the second term survives because ω_(α) _(i) _(β) _(i) =ω_(αβ) and there are N−1 terms. When the calculation of Eq. (30) is repeated for other transitions, one may find that ρ_(αβ) appears one time with coefficient (N+2)γ/2 and N−1 times with coefficient −γ/2. Therefore when the summation over the whole set is performed, one may find that the coefficient of ρ_(αβ) is (N+2)γ/2−(N−1)γ/2=3γ/2 which is exactly the linewidth of the independent atoms. However, if all dipole dipole interaction cannot be neglected, then the second term goes away in the secular approximation due to ω_(α) _(i) _(β) _(i) ≠ω_(αβ) and the linewidth in this case is (N+2)γ/2.

In Eq. (25), spontaneous emission terms (i≠j) have also been correlated. These terms appear only when the dipole-dipole interaction cannot be neglected. If an extreme case wherein all dipole dipole interaction cannot be neglected and no spectrum is overlapped, one has:

$\begin{matrix} \begin{matrix} {{S_{j}^{+}S_{i}^{-}\left. \alpha \right\rangle} = {{\frac{1}{\sqrt{2^{N}}}\underset{{k \neq i},j}{\otimes}\left( {\left. g_{k} \right\rangle + {\chi_{k}^{\alpha}\left. e_{k} \right\rangle}} \right)}\left( {\chi_{i}^{\alpha}\left. g_{i} \right\rangle} \right)\left. e_{j} \right\rangle}} \\ {= {{\frac{1}{4}\chi_{i}^{\alpha}\chi_{j}^{\alpha}\left. \alpha \right\rangle} + \ldots}} \end{matrix} & (29) \\ \begin{matrix} {{S_{j}^{+}S_{i}^{-}\left. \beta \right\rangle} = {{\frac{1}{\sqrt{2^{N}}}\underset{{k \neq i},j}{\otimes}\left( {\left. g_{k} \right\rangle + {\chi_{k}^{\beta}\left. e_{k} \right\rangle}} \right)}\left. e_{i} \right\rangle\left( {\chi_{i}^{\beta}\left. g_{i} \right\rangle} \right)}} \\ {= {{\frac{1}{4}\chi_{i}^{\beta}\chi_{j}^{\beta}\left. \beta \right\rangle} + \ldots}} \end{matrix} & (30) \\ \begin{matrix} \left. {{S_{j}^{+}\left. \alpha \right\rangle} = {{\frac{1}{\sqrt{2^{N}}}\underset{k \neq j}{\otimes}\left( {\left. g_{k} \right\rangle + {\chi_{k}^{\alpha}\left. e_{k} \right\rangle}} \right)}\left. e_{j} \right\rangle}} \right) \\ {= {{\frac{\chi_{i}^{\alpha}}{2}\left. \alpha \right\rangle} + \ldots}} \end{matrix} & (31) \\ \begin{matrix} {{S_{i}^{+}\left. \beta \right\rangle} = {{\frac{1}{\sqrt{2^{N}}}\underset{k \neq i}{\otimes}\left( {\left. g_{k} \right\rangle + {\chi_{k}^{\beta}\left. e_{k} \right\rangle}} \right)}\left. e_{j} \right\rangle}} \\ {= {{\frac{\chi_{i}^{\beta}}{2}\left. \beta \right\rangle} + \ldots}} \end{matrix} & (32) \end{matrix}$ where . . . denote terms that have different transition frequencies from ω_(αβ) and they can be neglected in the secular approximation. Then:

$\begin{matrix} {\sum\limits_{i \neq j}{\gamma_{ij}\left( {\quad\left\langle {A{{{S_{i}^{+}S_{j}^{-}\rho\left. \beta \right\rangle} + \left\langle {\alpha{{{\rho\; S_{i}^{+}S_{j}^{-}\left. \beta \right\rangle} - {2\left\langle {{\alpha\left. {S_{j}^{-}\rho\; S_{i}^{+}\left. \beta \right\rangle} \right)} = {{{\sum\limits_{i \neq j}{{\frac{{\gamma\;}_{ij}}{4`}\left\lbrack {{\chi_{i}^{\alpha}\chi_{j}^{\alpha}} + {\chi_{i}^{\beta}\chi_{j}^{\beta}} - {2\chi_{j}^{\alpha}\chi_{i}^{\beta}}} \right\rbrack}\rho_{\alpha\beta}}} + \ldots} = {{{\left\{ {{\sum\limits_{i \neq j \neq p}{\frac{\gamma_{ij}}{4}\left\lbrack {{\chi_{i}^{\alpha}\chi_{j}^{\alpha}} + {\chi_{i}^{\beta}\chi_{j}^{\beta}} - {2\chi_{j}^{\alpha}\chi_{i}^{\beta}}} \right\rbrack}} + {\sum\limits_{j \neq p}{\frac{\gamma_{pj}}{4}\left\lbrack {{\chi_{p}^{\alpha}\chi_{j}^{\alpha}} + {\chi_{p}^{\beta}\chi_{j}^{\beta}} - {2\chi_{j}^{\alpha}\chi_{p}^{\beta}}} \right\rbrack}} + {\sum\limits_{i \neq p}{\frac{\gamma_{ip}}{4}\left\lbrack {{\chi_{i}^{\alpha}\chi_{p}^{\alpha}} + {\chi_{i}^{\beta}\chi_{p}^{\beta}} - {2\chi_{p}^{\alpha}\chi_{i}^{\beta}}} \right\rbrack}}} \right\}\rho_{\alpha\beta}} + \ldots} = 0}}} \right.}}}} \right.}}} \right.} \right.}} & (33) \end{matrix}$

The first summation vanishes because χ_(i) ^(α)=χ_(i) ^(β) and χ_(j) ^(α)=χ_(j) ^(β) for i,j≠p. Because χ_(p) ^(α)=−χ_(p) ^(β) and χ_(j) ^(α)=χ_(j) ^(β), the first two terms in the second summation drop out. Similarly, the first two terms of the third summation also drop out. The remaining terms in the second summation and the third summation are just opposite because

$\begin{matrix} {{\sum\limits_{j \neq p}{\frac{\gamma_{pj}}{4}\left\lbrack {{- 2}\chi_{j}^{\alpha}\chi_{p}^{\beta}} \right\rbrack}} = {{\sum\limits_{i \neq p}{\frac{\gamma_{pi}}{4}\left\lbrack {{- 2}\chi_{i}^{\alpha}\chi_{p}^{\beta}} \right\rbrack}} = {\sum\limits_{i \neq p}{\frac{\gamma_{pi}}{4}\left\lbrack {{+ 2}\chi_{i}^{\beta}\chi_{p}^{\alpha}} \right\rbrack}}}} & (34) \end{matrix}$

When the first order correction of the wavefunction is considered, this correlated spontaneous emission rate is nonzero. However their values are on the order of (Ω_(ij)/Ω_(i))² which may be assumed to be very small.

Although particular embodiments of the present invention have been shown and described, they are not intended to limit what this patent covers. One skilled in the art will understand that various changes and modifications may be made without departing from the scope of the present invention as literally and equivalently covered by the following claims. 

What is claimed is:
 1. A method for determining the positions of a plurality of two-level atoms located in a line comprising: shining two, linear-polarized laser fields on atoms trapped in a quantum dot from opposite directions such that a standing wave is formed; monitoring the resonance-fluorescence photons emitted by the trapped atoms; obtaining the Rabi frequencies from the sideband peaks in the resulting fluorescence spectrum; calculating the positions of the trapped atoms from the Rabi frequencies.
 2. A method as recited in claim 1 wherein the atoms are identical atoms.
 3. A method as recited in claim 1 wherein the atoms are located within one wavelength of the laser fields.
 4. A method as recited in claim 1 further comprising stretching the standing wave such that the sample is located within the linear region between the node and antinode.
 5. A method as recited in claim 1 further comprising detecting fluorescence photons with a detector in the z direction when the atoms are located in a line along the x axis.
 6. A method as recited in claim 1 wherein far-field optical microscopy is used to obtain the resonance fluorescence spectrum.
 7. A method for determining the positions of a number of two-level atoms in two dimensions (x,y) comprising: shining a first gradient laser field along the x-axis direction; obtaining a discrete set of x-axis position information of the atoms from the resulting first resonance fluorescence spectrum; shining a second gradient laser field along the y-axis direction; obtaining a discrete set of y-axis position information of the atoms from the resulting second resonance fluorescence spectrum; shining a third gradient field from a direction which avoids any two pairs of (x,y); and, determining the positions of the atoms from the resulting third resonance fluorescence spectrum.
 8. A method as recited in claim 7 wherein far-field optical microscopy is used to obtain the resonance fluorescence spectra.
 9. A method for determining the positions of a plurality of two-level atoms located in a line comprising: shining two, linear-polarized laser fields on atoms trapped in nitrogen vacancy centers in diamond from opposite directions such that a standing wave is formed; monitoring the resonance-fluorescence photons emitted by the trapped atoms; obtaining the Rabi frequencies from the sideband peaks in the resulting fluorescence spectrum; calculating the positions of the trapped atoms from the Rabi frequencies.
 10. A method for determining the positions of a plurality of two-level atoms located in a line comprising: shining two, linear-polarized, gradient laser fields on the atoms from opposite directions such that a standing wave is formed; monitoring the resonance-fluorescence photons emitted by the atoms; obtaining the Rabi frequencies from the sideband peaks in the resulting fluorescence spectrum; calculating the positions of the atoms from the Rabi frequencies.
 11. A method as recited in claim 10 further comprising changing the gradient of the laser field by a certain amount.
 12. A method as recited in claim 10 further comprising applying a first, weaker gradient field and subsequently applying a second, stronger gradient field.
 13. A method as recited in claim 12 wherein the second, stronger gradient field is less than the dipole-dipole interaction energy of the atoms.
 14. A method for determining the positions of a plurality of two-level atoms located in a line comprising: shining two, linear-polarized laser fields on the atoms from opposite directions such that a standing wave is formed; monitoring the resonance-fluorescence photons emitted by the atoms; obtaining the Rabi frequencies from the sideband peaks in the resulting fluorescence spectrum; calculating the positions of the atoms from the Rabi frequencies; and, averaging the results when two results are obtained for the positions of each atom. 